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Abstract. Recent RHIC results on pion production in AuAu collision at y/s = 130 and 200 AGeV display 
a strong suppression effect at high pr- This suppression can be connected to final state effects, namely 
jet energy loss induced by the produced dense colored matter. Applying our pQCD-based parton model 
we perform a quantitative analysis of the measured suppression pattern and determine the opacity of the 
produced deconfined matter. 

PACS. 12.38.Mh Quark-gluon plasma in quantum chromodynamics - 24.85+p Quarks, gluons, and QCD 
in nuclei and nuclear processes - 25.75-q Relativistic heavy-ion collisions 



1 Introduction 

The experimental data on high-p^ tt° production in cen- 
tral AuAu collisions at mid-rapidity at y/s — 130 and 200 
AGeV have shown a strong suppression compared to bi- 
nary scaled pp data ■ This suppression vanishes with 
increasing centrality and no effect appears in peripheral 
AuAu collisions. Thus a detailed quantitative analysis of 
the suppression pattern ( "jet tomography" jjj) yields infor- 
mation about the properties of the produced dense mat- 
ter and the impact parameter dependence of the forma- 
tion of deconfined matter in AuAu collisions. Recent data 
on dAu collisions [4"||^] validate our effort: since no sup- 
pression was found at mid-rapidity in the dAu reaction, 
initial state effects (e.g. a strong modification in the in- 
ternal parton structure of the accelerated Au nuclei) can 
not be responsible for the measured suppression in AuAu 
collisions. Thus, induced jet-energy loss in the final state 
becomes a strong candidate to explain the missing pion 
yield. Jet energy loss can be calculated in a perturbative 
quantum chromodynamics (pQCD) frame [HUZj. 

We investigate jet energy loss in a pQCD improved 
parton model. In Sect. |2] we introduce the basis of our 
model, especially a phenomenological intrinsic transverse 
momentum distribution (intrinsic kx) for the colliding nu- 
cleons, which is necessary to reach a better agreement be- 
tween data and calculations in pp collisions 8 9 JlOllllj . 
For nucleus-nucleus (AA) collisions, initial state effects are 
considered, e.g. nuclear multiple scattering, saturation in 
the number of semihard collisions, and a weak shadowing 
effect inside the nucleus |51HTj) . In Sect.|3|we summarize 
the description of induced gluon radiation in thin non- 



Abelian matter and include the GLV-description [J] of en- 
ergy loss into the pQCD improved parton model. This way 
our model becomes capable to extract the opacity values 
of the produced colored matter at different centralities. In 
Sect. 0] we discuss the obtained results. 



2 Initial State Effects in AuAu Collisions 

The invariant cross section of pion production in an AA' 
collision can be described in a pQCD-improved parton 
model developed for pp collision and extended by a Glauber- 
type collision geometry and initial state nuclear effects for 
A A' collisions as [T31ll4lll2| : 
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Here £^(6) = J dzpA(b,z) is the nuclear thickness func- 
tion normalized as J d 2 btA(b) = A. For small nuclei we 
use a sharp sphere approximation, while for larger nuclei 
the Wood-Saxon formula is applied. 
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In our next-to-leading order (NLO) calculation 12 , 
da/dv represents the Born cross section of the partonic 
subprocess and K a b iC (s, v, w, Q, Q r , Q) is the correspond- 
ing higher order correction term, see Ref.s |12U13U14| . We 
fix the factorization and renormalization scales and con- 
nect them to the momentum of the intermediate jet, Q — 
Q r — (4/3)p q (where p q = pt/z c ), reproducing pp data 
with high precision at high pt ^Hl ■ 

The approximate form of the 3-dimcnsional parton dis- 
tribution function (PDF) is the following: 

fa/p(x a , k Ta ,Q 2 ) = fa/p(x a ,Q 2 ) ■ ga/pi^Ta) ■ (2) 

Here, the function f a / P (x ai Q 2 ) represents the standard 
longitudinal NLO PDF as a function of momentum frac- 
tion of the incoming parton, x a at scale Q (in the present 
calculations we use the MRST(cg)[in] PDFs). The par- 
tonic transverse-momentum distribution in 2 dimensions, 
9a/pO^T), is characterized by an "intrinsic fcr" parameter 
as in Ref.s |10U12| . In our phenomenological approach this 
component is described by a Gaussian function |1()U11| . 

Nuclear multiscattering is accounted for through broad- 
ening of the incoming parton's transverse momentum dis- 
tribution function, namely an increase in the width of the 
Gaussian: 



C ■ h pA (b) 



(3) 



Here, (k^)p P = 2.5 GeV 2 is the width of the transverse 
momentum distribution of partons in pp collisions |1()U15| , 
h p A(b) describes the number of effective NN collisions at 
impact parameter b, which impart an average transverse 
momentum squared C. The effectivity function h p A(b) can 
be written in terms of the number of collisions suffered by 
the incoming proton in the target nucleus. In Ref. ^Oj 
we have found a limited number of semihard collisions, 
v m = 4 and the value C — 0.4 GeV 2 . 

We take into account the isospin asymmetry by us- 
ing a linear combination of p and n PDFs. The applied 
PDFs are also modified inside nuclei by the "shadowing" 
effect [HI- 

The last term in the convolution of eq. Q is the frag- 
mentation function (FF), D™(z c , Q 2 )- This gives the prob- 
ability for parton c to fragment into a pion with momen- 
tum fraction z c at fragmentation scale Q = (4/3)pt- We 
apply the KKP parametrization |18| . 
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Fig. 1. Pion production in central AuAu collision with the 
calculated opacities n = 3.5 ± 0.25 (upper panel). Data are 
from PHENIX Collaboration |T]. The lower panel displays a 
comparison between the data and the calculations. 



given by the following form |7|: 
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where Cr is the color Casimir of the jet, fJ./X g ~ a 2 p pa rt 
is a transport coefficient of the medium, proportional to 
the parton density, p par t- The color Debye screening scale 
is /i, and X g is the radiated gluon mean free path. N(E) 
is an energy dependent factor with asymptotic value 4. 

Considering a time-averaged, static plasma, the aver- 
age energy loss, AE, will modify the argument of the FFs: 



3 Jet-Quenching as a Final State Effect 

The energy loss of high-energy quark and gluon jets trav- 
eling through dense colored matter is able to give informa- 
tion on the density of gluons jHJ- This non-Abelian radia- 
tive energy loss AE(E, L) can be described as a function 
of gluon density: h = L/X gi the mean number of jet scat- 
terings, where L is the length traversed by the jet and 
X g is the mean free path in non-Abelian dense matter. In 
"thin plasma" approximation energy loss in first order is 
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Here z* = z c / (1 — AE/p c ) is the modified momentum 
fraction. 

In Fig. 1 we present our result on pion production 
in most central (0 — 10%) AuAu collisions at yfs = 200 
AGeV. We included all initial and final state effects dis- 
cussed above and used the opacity value n = 3.5 ± 0.25 to 
reproduce the experimental data. 
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Fig. 2. Upper panel displays the pion production in AuAu col- 
lision in different centrality bins with the calculated opacities. 
Data are from PHENIX Collaboration 1 . Lower panel shows 
the pion production in dAu collision |1| and the result of our 
calculation |f 5| . 



4 Centrality Dependence in AuAu collisions 

In the top panel of Fig. [3 we display the nuclear modifi- 
cation factor 



Raa'(pt, b) = 



1 E n da^ A (&)/dV 
N bm ' EM P /d 3 pr 



(0) 



(where Nun is the number of binary collisions), as a func- 
tion of pt at different impact parameter ranges. The mea- 
sured suppression is reproduced in the most central col- 
lisions with opacity n = 3.5 ± 0.25 as it was shown in 
Fig-ffl Curves with lower values of opacities are shown for 
comparison to the more peripheral data. Since the central- 
ity 60 — 92% includes non-peripheral events, the obtained 
n = 1.5 seems to be a reasonable opacity value for this bin. 
In the very peripheral case the opacity should be reduced 
to a small value ^H], however ±20% error on the nuclear 
modification factor does not allow more detailed investi- 
gation. A more quantitative analysis can be performed, 
when data will be available with smaller error bars. At 
higher precision, the geometry of the hot overlap zone can 
be included and a more detailed analysis of the impact 
parameter dependence can be accomplished. In this case 
the properties of the produced hot matter will be studied 
in a more quantitative way. 

In the bottom panel of Fig. [21 the dAu data are com- 
pared to our calculation JS] to demonstrate the validity 
of our description for the initial nuclear effects. 
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